Quantum- Classical Correspondence via Liouville Dynamics: II. 
Correspondence for Chaotic Hamiltonian Systems 



Joshua Wilkie and Paul Brumer 
Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, 
)P ' Ontario, Canada M5S lAl 

< 

in 



Abstract 



We prove quantum-classical correspondence for bound conservative classi- 

^ ■ cally chaotic Hamiltonian systems. In particular, quantum Liouville spectral 

■ projection operators and spectral densities, and hence classical dynamics, are 

00 

(<~^ ■ shown to approach their classical analogs in the /i ^ limit. Correspondence 
^ ■ 

0^ . is shown to occur via the elimination of essential singularities. In addition, ap- 

J>^. plications to matrix elements of observables in chaotic systems are discussed. 

o 

^ : I. INTRODUCTION 

o 

> 

\ The validity of quantum mechanics as a description of the macroscopic world is contin- 

\ gent upon the reduction of the laws of quantum mechanics to Newton's laws in the limit 
where the characteristic actions of a system are large with respect to Planck's constant [|I|. 
Thus, diagonal and off-diagonal matrix elements must reduce to their classical analogs and 
quantum dynamics must reproduce the predictions of classical mechanics as /i — > 0. Despite 
the fundamental importance of quantum-classical correspondence it has only been satisfac- 
torily demonstrated 0- in the very restrictive case of regular systems, i.e., systems which 
classically possess as many constants of the motion as degrees of freedom. Indeed some au- 
thors have suggested that bound quantum systems with discrete quantum spectrum and 
chaotic classical analog may violate the correspondence principle 0. These doubts about 
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the validity of the correspondence principle for chaotic systems stem from the difficulty of 
reconciling the quasiperiodic nature of bound state quantum dynamics with the chaotic 
nature of classical dynamics for the same Hamiltonian. The issue of correspondence for 
quantum systems whose classical analogs exhibit chaos (irregular systems) is thus of great 
interest . 

Verification of correspondence should be distinguished from the development of semiclas- 
sical approximation methods. While semiclassical theories provide a natural starting point 
for an exploration of the classical limit their existence does not guarantee correspondence. 
For example, semiclassical theories for regular systems preceded the development of modern 
quantum mechanics , but an understanding of correspondence for regular systems has only 
recently been achieved By comparison, attempts to develop semiclassical quantiza- 

tion rules for chaotic systems have had some success |^, whereas the correspondence limit 
remains largely unexplored . In this paper we demonstrate that the existing semiclassical 
theories of quantum dynamics for classically chaotic systems are sufficiently well developed 
to allow us to show that such systems do in fact approach their proper correspondence limits 
as Planck's constant approaches zero. This completes the Liouville correspondence program 
outlined in the preceding paper and significantly extends the results of our study of 
quantum maps where we rigorously demonstrated that a nonchaotic quantum map 

dynamics can completely recover a fully chaotic classical dynamics in the limit h —>■ 0. 

The Liouville picture affords a means of gaining insight into the connections between 
quantum and classical mechanics and is a natural framework for studies of corre- 

spondence. As outlined in the preceding paper (henceforth referred to as Paper 1) the 
essential ingredients for Liouville dynamics are eigenstates and eigenvalues of the Liouville 
operators in both mechanics. In particular, the dynamics is completely characterized by the 
Liouville eigenfunctions and eigenvalues or the spectral projectors once the class of allowed 
initial distributions is specified. Here we consider correspondence in chaotic systems from 
this Liouville perspective. 

Quantum Liouville eigenfunctions for conservative Hamiltonian systems whose classical 
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analogs are chaotic take the form |n)(m| where \n) are eigenstates of the Hamiltonian, i.e., 
H\n) = En\n). These distributions are eigenfunctions of the complete set of operators I/,7tI 
where L = \[H, ] is the quantum Liouville operator and where Ti. = ^[H, ]+ is the Hermitian 
energy operator in the Liouville picture 0. That is, they are solutions of both the time 
independent Liouville equation 

L\n){m\ = \n,m\n){m\, (1) 
where Xn,m = {En — Em)/h, and of the energy eigenequation 

H\n){m\ = En,m\n){m\, (2) 

with En,m = {En + Em)/ 2. 

Consistent with von Neumann's criteria for quantum ergodicity [0 , we deal with quan- 



tum systems with chaotic classical analog |15jT6[| for which the spectrum of energies E, 



is nondegenerate. For such systems the states \n){m\ are specified by the integers n and 
m, or equivalently by the frequency \n,m and energy En^m- Since the distributions \n){m\ 
govern the quantum dynamics ^ an understanding of their h —>■ limit, or of their Wigner 
representation p^^(x), 

PlJ^) ^ h'^^' J rfve^P-^/^(q-v/2|n)(m|q + v/2), (3) 

would seem essential for verification of correspondence. [Here x = (p, q) where p are 
the momenta and q are the coordinates.] However, as shown below, the relevant objects 
for the study of correspondence in chaotic systems are the quantum spectral projection 
operators '^7^ which are of the form p5^^(xo)p^.^(x) in the Wigner representation. That 
is, we demonstrate that for irregular systems these quantum Liouville spectral projection 
operators approach classical spectral projection operators T of the same frequency and 
energy as /i — 0, i.e., that 

pZi^o)pln{^) - dE T^„(x,xo), (4) 

and 



Pn!mi^o)pn,mi^) ^ dEdX Ti5„.„,A„,™ (x; Xq), {u ^ Til) . (5) 

Here the distributions T^; and T^; ;^ are the stationary and nonstationary Liouville spectral 
projection operators of classical dynamics ||^, discussed in Paper 1. We also show that the 
spectrum of the quantum Liouville operator goes to that of the classical operator as h —>■ 
and that the correspondence emerges smoothly via the elimination of essential singularities. 
Proof of Eqs. and (^), plus proof of the correspondence of the Liouville spectra, suffices 
to prove quantum-classical correspondence in chaotic systems. 

Note that, unlike their quantum analogs p^!„(xo)p^„(x), the nonstationary chaotic clas- 
sical spectral projection operators '^E„,„,Xnm cannot be written as a product of Liouville 
eigenfunctions but rather consist of a sum of products of Liouville eigenfunctions (see Paper 
1). Hence Eq. suggests, as discussed below, that the individual nonstationary quantum 
Liouville eigenfunctions p^^(x) for quantum systems with chaotic classical analogs do not 
have well defined correspondence limits. This situation is quite different from that of the 
integrable case discussed in Paper 1 and necessitates the introduction of new tools to prove 
correspondence. 

This paper is organized as follows: Section ^ introduces a useful Dirac notation to 
simplify our formal manipulations, and the proof of Eqs. and is expressed in this 
Dirac form. Section |T| proves correspondence for both the Liouville spectral projectors 
and the Liouville eigenvalues. This treatment ignores higher order corrections relating to 



scars, which are treated in Section The proof of correspondence allows us to consider 
the classical limit of operator matrix elements, which is discussed in Section 0. Section ^ 
provides a summary. 



II. A DIRAC FORMULATION OF LIOUVILLE DYNAMICS 

The effectiveness of the Liouville picture is limited by the clumsiness of the associated 
density matrix notation. In this section we introduce a useful Dirac notation which simpli- 
fies manipulations considerably We will also employ a Dirac notation for the classical 
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Liouville dynamics in order to maintain symmetry between the quantum and classical for- 
mulations. 

Let \n) be a complete, orthonormal set of basis states for the quantum Hilbert space 
associated with the solutions of the Schrodinger equation. That is, {n\m) = 5n,m and 
Y.n = 1- From these states we construct distributions Pn,m = \n){'m\ which are a basis 

in the Hilbert space associated with solutions of the von Neumann equation. It is natural 
to assign a Dirac notation to these basis states, i.e., 

\n,m) = Pn,m- (6) 

A complete orthonormal basis \n) of Schrodinger states then yields a complete set of Liouville 
states \n,m). One can now easily deduce that the dual space is spanned by the linear 
functionals 

{n,m\=TT{pi„^-} (7) 

by requiring that {n,m\k,l) = Sn,m^k,i- Note that the normalization of the states \n,m) has 
been chosen so that \n,m){n,m\ is a projection operator. Completeness implies that 

\n,m){n,m\ = 1. (8) 

n,m 

The spectral decomposition of the Liouville operator then takes the form 

L = ^Xn,m\n,m){n,m\. (9) 



The Hermitian operators \n,m){n,m\ = p„ „iTr{/)Jj ^-j are obviously the spectral 



projection operators in the Liouville picture, in the same way that \n){n\ are the spectral 
projection operators in the traditional Hamiltonian picture. Arbitrary superoperators of the 
form [O, ]± (i.e., [O, ]±p = Op ± pO), of interest below, can be expanded on the \n,m) 
states as 

[0,]±= E \n,m){n,m\[d, U\k,l){k,l\ (10) 

n,m,k,l 



where the superoperator "matrix elements" are (n, m|[0, ]±|A;,/) = On^k^i^m ± Oi^rn^n,k- 
Physical states \p) are defined as 

\p) = h-^/'p, (11) 

with corresponding kets 

(p| = h-^/'Tr{p^.} = h-^/'TT{p.}, (12) 

with the latter equality due to the fact that p'^ = p. Equations ([TT|) and (|12D, which define 
the physical states, differ from Eqs. (|]) and (^ which define the basis states by a factor 
of /i"''/^ which is introduced so that the quasi-probability distributions associated with the 
physical states have the correct dimensions in a phase space representation, i.e., inverse 
action to a power equal to the number of degrees of freedom. We also assign states \A) to 
operators A (i.e., operators operating on the Hilbert space spanned by \n)) via \A) = h^/'^A 
and {A\ = h''/'^TT{A'^ ■}. The expectation of A is then given by {p\A) = {A\p)* = Tt{pA}. 
In this notation the von Neumann (quantum Liouville) equation is 

||p(t)) = ~^L\p{t)). (13) 



and the Wigner-Weyl representation [0,0 of a state p takes the form 



(x|p) = h-'/^TT{h-'/^A{^)p}, (14) 

where 

A(x) = h-' J dudw e*[^-(p-p)+"-(i-q)]/^. (15) 



Thus, employing Eq. (PA|) , we identify 

(x| = Tr{/i-^/2^(x)-} = h-'/^Tr{A\^)-} (16) 

where the second equality is due to the fact that A(x) is Hermitian. The particular form 
of the corresponding ket [ P5| , P^ is determined by demanding that (x|x') = (5(x — x'). Thus, 
here 
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|x) = /i-^/2A(x), (17) 

where (x|x') = /t~*Tr{A(x)A(x')} = 5(x — x'). Since |x) and (x| span the Hilbert space 
and its dual space, they satisfy the closure relation 

J rfx |x)(x| = 1. (18) 

Definitions (|TBp and (|T^ in conjunction with Eqs. (|TT]) and (|T2|) guarantee that the prob- 
ability densities (x|p) have the correct dimensions. Other phase space representations [in 
which (x| and |x) may be quite dissimilar], and a general transformation theory between 
them is provided elsewhere |^ . 



Consider then the Liouville spectral decomposition [i.e., Eqs. (|T]) and (|^)] for a chaotic 
quantum system in the Dirac notation. As eigenfunctions of L and Ti, the \n,m) satisfy 

L\n,m) = Xn,m\n,m), (19) 

and 

H\n,m) = Enrn\n,m). (20) 
In the Wigner-Weyl representation Eqs. (|19D and (|20|) become 

(x|I/|ra,m) = L(x)(x|n,m) = A„,m(x|n, m), (21) 

and 

(x|7i|ra,m) = 7i(x)(x|r2, m) = £'„_m(x|n, m), (22) 
where ^^(x) = ^if(x) sin(/icr/2) is the quantum Liouville operator, and 7Y(x) = 



n 

H('x.) cos{ha/2} is the energy operator. Here = ^ ■ J ■ ^ is the Poisson bracket, i.e.. 



A{x)aB{x) = {A,B}, and J 



lo-A 



yl J 



is the 2s x 2s dimensional symplectic matrix 



27[ . Expanding ?i(x) in powers of h shows that the classical analog of 7^(x) is the energy 



function H{x.), and that Lc(x) is the correspondence limit of -^v(x). 



Consider now the classical case. The classical analog of a phase space representation is a 
choice of canonical variables for a classical distribution p^- Thus we denote the phase space 
representation of pc by Pc(x) = (x|pc)- 

The classical Liouville spectral decomposition, and the properties of the eigendistribu- 
tions discussed in Paper 1 |^ are readily restated using the Dirac notation. Associating 
states \E) with the classical distributions Pe{^) which span the point spectrum, and states 
\E,X,i) with the classical distributions p|;;s^(x) which span the continuous spectrum, the 
full set of equations for the spectral decomposition becomes: 

{E'\E)=5{E'-E), (23) 
{E'\E,\,i) = 0, (24) 
(E'X', i'\E, A, i) = 6e,eS{E' - E)6{X' - A), (25) 

roo roo r 

dE\E){E\+J^ dEjdXj2\E,\,i){E,X,i\ = l, (26) 

e-'^^^'lE) = \E), (27) 

and 

e-'^^'\E,X,i) =e-'^'\E,X,i). (28) 

Here the line through the integral in Eq. (p6|) indicates that the point spectrum eigenvalue 
A = has been removed (see Paper 1). Two further equations relate to the second constant 
of the motion, a classical energy operator He- 

n,\E) = E\E), (29) 

and 

HclE^XJ) = E\E,X,e). (30) 
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In the phase space representation parameterized by x these equations become 

(x|K|^) = H{^){^\E) = E{^\E), (31) 

and 

(x|K|^, A, i) = i7(x)(x|E, A, i) = E{x\E, A, i). (32) 

A complete set of stationary and nonstationary classical Liouville eigenf unctions p_e;(x) = 
{'x.\E) and p^;^(x) = {yi\E,\,£), were introduced in Paper 1 where the integer i labels the 
infinite degeneracy of the continuous spectrum [|l^. In addition, spectral projection opera- 
tors T£;(x;xo) and Te,\{'X',^o) were introduced; these are the phase space representations 
of the classical operators 6 {E -He) = \E){E\ and 6{E -Hc)S{X- L^) = Z^lE, X,i){E, X,i\. 
Specifically, 

(x|E)(E|xo) = Te{^] xo) = Jdx'6{E - H{:>c')) " PeWPe{^o) (33) 

and 

(x|E,A,£)(E,A,£|xo)=T£;,;,(x;Xo) 

e 

= -6{E - if(xo)) / dt'e^'' <5(xo - X(x, -t')) = ^ p%,{^o)p'E,xi^), (34) 

where X(x, —t') is the phase space point from which x emerges over a time t' . 

In terms of these eigenfunctions the spectral decomposition [Eq. (|26|) ] takes the form 

/•oo roc r 

dE p%{^o)pe{^) + dEjdX Ept,A(xo)p'E,A(x) = 5(x - xo), (35) 

or 

roo roo r 

I c/E Ti5(x; Xo) + / c/E-fc/A T£;a(x;xo) = 5(x-xo). (36) 

JO JO J 

Thus, the evolution of any initial distribution p(x, 0) can be written, in quantum me- 
chanics, as an expansion: 
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I (ixop(xo, 0) E p-;(xo)p-„(x) + E p::Mplmi^)e-''--'] (37) 



and in classical mechanics as 



P 



(x,t) = J dE CE Pe{^) + J dE jdXj2cE,x,i pE,x{^)e''^' 

= J d^op{^o,0)[J dE pU^o)pe{^) + J dE-fdXj2pE,xMPEA^)e-'^'] 

= J rfxop(xo, 0)[JdE Ts(x; xq) + J dE -fdXTE,x{x; xo)e-'^*]. (38) 



Equations (|37| ) and ( pSf ) make clear that a demonstration of correspondence for the spectral 
projection operators and their eigenvalues is sufficient to establish correspondence for the dy- 
namics, i.e., that p(x, t) —>■ p'^(x, t) as h —>■ 0. That is, formally establishing correspondence 
requires demonstrating 

|n,n)(n,n| -> dE \En)iEn\ (39) 

and 

\n,Tn){n,m\ ^ dEdX J2 \En,m, Xn,mJ)iEn,m, Xn,mJ\, (40) 

i 

or 

(x|n, n)(n, n|xo) = p:,;(xo)p:^, Jx) ^ T^Jx; xq) = (x|E„)(Ejxo) (41) 

and 

(x|n,m)(n,m|xo) = P^,tn(xo)p;^,^(x) ^ dEdX T^^ ^,a„,,„(x; xq) 

= dEdX E (x|-E„^m, A„,^m, £)(-E„_m, A„_m, £|xo), (42) 



with the infinitesimals dE and dX to be determined. These limits are proven in Sec. [Ill A 
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III. CORRESPONDENCE 



Consider then the correspondence hmit, i.e., the hmit of the quantum Liouville dynamics 
as h 0, with the h ^ hmit taken before the T oo hmit |^,28|. This order, h ^ 



first, is consistent with the actual physics in which one first chooses a particular system 
and then propagates it for long times. Technically, this is achieved by first broadening the 
system energy by some amount e ^ h/T^i^ (thus restricting the dynamics to finite time) 
and taking the h ^ limit with e fixed. The broadening e can then be chosen infinitesimal, 
e — s> (AE), (where (AE) is the average spacing between neighboring energy levels) allowing 
for long time dynamics. Here Tmin is the period of the shortest periodic orbit. 

The physical significance of correspondence under these limits is clear. A transition from 
quantum to classical behavior will be observed in the dynamics of a physical system a.s h —>■ 
provided that (a) the apparatus with which we observe its dynamics has a fixed, classically 
small but quantum mechanically large, energy resolution, and that (b) we do not observe 
its dynamics beyond the recurrence time given approximately hj h/ (AE). 

A. Correspondence for Spectral Projection Operators 

Here we examine the correspondence limits of the spectral projection operators 
\n,n){n,n\ and \n,m){n,m\. We focus attention on the nonstationary case. The stationary 
case [Eq. (^Tj)] has already been obtained by Berry and Voros but we work through 

this case to demonstrate the consistency of our approach. In the latter case consider Berry's 
formula iiO,|3ll 



lini7reiy(x;E„,e) = J dq'e'^-^^'^^iq- q' /2\n){n\q + q' /2) = h'''^ (x|n,n) (43) 

for finite e, in order to investigate the semiclassical form of the stationary Liouville eigendis- 
tribution (x|n, n) in the Wigner representation. Here (x; E, e) is the Lorentzian weighted 
sum of Wigner functions over a width e about an energy E, i.e.. 
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W{^-E,e) = J ciq'e-*P-'i'/^(q + q72|(5,(E-^)|q-q72), (44) 



6,{E -H) = --Im \ = ^ j ' dt e'^E-H)tin^-e\t\in_ ^45^ 



where 

sAF.-m = - 

vr E-H + ie h 

For quantum systems with chaotic classical analogs, Berry ||3^ has shown, for small h and 
small e, that (where ~ denotes the form in the limit) 

W^(x; e) ~ 6,iE - H{^)) + ^ iy^,,,(x; E, e), (46) 

j 

and thus, by employing Eq. (^31), that 

(x|n, n) ~ neh-'/\iEr, - //(x)) + ireh-'/^ ^ iyi,,,(x; e). (47) 

j 

Here M/^car is of the order h^~^ smaller than the 6e{En — H{x.)) term, and hence vanishes 



rapidly as h ^ 0. These scar terms, neglected in this section, are considered in Sec. |y. 
Neglecting the scar terms gives 

(x|n, n) ~ h-'/\e6{En - H{^)) ~ 7re{x\En) / {AE^/^ (48) 

where (x|E„) = ps„(x) = 5(E„ - /J(x))/[/ dx'5(E„ - H{^'))Y/^ (see Paper 1) and where 
(AE) = / J d'x'5{E — H(x.')) is the average adjacent energy level spacing. Therefore 

(x|n,n)(n,n|xo) ~ (7re)2(x|^„)(K|xo)/(A^). (49) 

Correspondence for the stationary eigenstates [Eq. (|4l|) ] then results if we take the limit 
vre (AE) and note that (A^) c/E. 

The proof of Eq. (|4^) follows in a similar fashion from the following important rela- 
tionship (proven in Appendix A) between stationary and nonstationary quantum Liouville 
eigenfunctions [p^] : 

(x|n,m)(n,m|xo) = h~' J rfuc/ve-^^P-P^^'^/^e^^'i-^")-"/^ • 

((p + u/2, q + v/2) |n, n) (m, m| (po - u/2, qo - v/2)) 
= h-' J rfye^(™)--^ y/'^ . (x + y/2\n, n) (m, m|xo - y/2). (50) 
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Here the tilde denotes the transpose (a column vector) of the row vector y. Substituting 
Eq. (H) into Eq. (0) gives: 

(x|n,m)(n,m|xo) = lim /i^^V^eiea / dy e^(''-^o)-J-y/?i . 

iy(x + y/2; En, ei)W*{xo - y/2; E^, e^) (51) 

if the energy eigenvalues En and Em are both nondegenerate, the case for a chaotic system. 
In the limit of small ei, €2 we obtain 

(x|n,m)(n,m|xo) ~ h-^'nhie2 J dy e*^™)-^-^/^ ■ 

+ y/2; ei)iy*(xo - y/2; es) (52) 

which is amenable to semiclassical analysis. 

As a first approximation we neglect the scar corrections and employ 

/oo 
dt e'iE~^i^))t/h^-e\t\/h (53^ 
-00 

as we did in the stationary case. Substituting Eq. (^) into Eq. ([5^), followed by a simple 
change of variables y ^ hy, yields 

(x|n,m)(n,m|xo) ~ Tr^eiCa / dy e=^"^("-"°)-^-^5e,(K, - H{^ + hy/2))6,,{Em - H{^o - hy/2)) 
J 

exp{i(K - + hy/2))ti/h} exp{2(E„ - i/(xo - hy/2))t2/n}. (54) 

Note, at this stage, the presence of essential singularities in each of the highly oscillatory 
phase factors ex]i{i{En — H{yi + hy /2))ti/Ti} and exp{i(£'m — H{-k.q — hy /2))t2/h}. 

We now let h with ei, €2 fixed. Expanding the displaced Hamiltonian functions in 
Eq. ( p^ in powers of h: 

H(. + hy/2) ^ ff(x) + 1^ . y ^ H(.) - J] -J-y (56) 

- hy/2) ~ m^) - . y ~ ff(x„) + J] 'J-J (56) 
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and substituting these expressions into Eq. ( p4D gives, 

exp{27rz ■ [(x - ^^Jt,/2) - (xo - jt,/2)] . J . y}. (57) 

The factor e~'^il*il/^e~'^2'*2'/^ guarantees that the integrand is zero for all but short times 
since ei\ti\/h ^ 27r|ti|/Tmm and e2\t2\/h ^ 27r|t2|/7'min, so that we can use the short time 



approximation: 



X(x, -ti/2) ~ X - Jti/2, (58) 

X'(xo, -h/2) ~ Xo - ^^^^Jt,/2. (59) 

9X0 



Using these results in Eq. gives 

(x|n,m)(n,m|xo) ~ / (iy(itidt2e-^i'*^'/"e-^^'*^'/"e'(^"-^W)ii A 

gi(i?„-/f(xo))t2/?ig^p|2^^[X(x, -ti/2) - X'(xo, -t2/2)] ■ J ■ y} 

gi(i?„-H(x))tiAgi(E^-j/(xo))t2/a^(X(x, -ti/2) - X'(xo, -t2/2)). (60) 

Noting that Eq. (^) is identically zero unless x and Xo are on the same trajectory and using 
the fact that the Hamiltonian is time independent allows us to replace -ff(x) by i/(xo) in 
the exponential. Next we perform a canonical transformation (i.e., time translation) to put 
Eq. (|60D in the form 

2 

(x|n,m)(n,m|xo) ~ I dtidt2e-'^\''\'^e-'^\'^\'^ ■ 

g*(i?„-H(xo))ti/7igi(£™-/f(xo))i2A^(x^ _ x(x, -(ti - t2)/2)) 



/l2 

^iiEn-H{xo))to/h^i{E^-H{xo)){to~2t' 



)/^5(xo-X(x,-t')). (61) 



where we have changed variables to t' = (ti — t2)/2 and to = ^i- Note that e~^^'*°'/^ ~ 

unless to ~ so that we can replace e-^2\to-2t'\/h -j^^ ^-2e2\t'\/h (pT]). Defining eo = ei/2 
and e' = 2e2 we obtain 
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(x|n,m)(n,m|xo) ~ / c^t'dtoe-'^°l*°l/^e-^'l*'l/^ 

^iiE„-H{yio))to/h^i{E„,-H(xo)){to~2t')/h^^^^ _ ^^^^ -t')). (62) 

Using identity (^5|) in reverse then yields two equivalent forms: 

2 / 

(x|n,m)(n,m|xo) ~ ^5.„(i?„,„. - i/(xo)) / cit'e-^'l*'l/^e-2^(^-"^(^«))*'/^5(xo - X(x, -t')) 

(63) 

and 

2 / 

(x|n,m)(n,m|xo) ~ - if(xo)) j dt'e-^'l*'l/^e^^"-*'(5(xo - X(x, -t')). (64) 

Note that Eq. (^) no longer exhibits the essential singularities present in Eq. (^4|). This 



is due to the expansions in Eqs. (|55D and (^) through which the essential singularities are 



eliminated. We found that the same mechanism, i.e., elimination of essential singularities. 



was responsible for correspondence in chaotic mappings of the torus [TQ • 

We now take the limit "T — * oo", that is we let eo/h,e'/h 0. Note that eo and e' 
essentially define a cutoff in time beyond which the semiclassical approximations break down. 
The commonly adopted cutoff time is the density of states time T^s ^ h/ (AE). Since our t' 
is symmetric about zero, propagation to T^s implies that — Tds/2 < t' < Tds/2. To achieve 
this we let e'T^s/^h — > 1 or, substituting Tds ~ h/{AE), vre' — > (AE). Thus the precise 
limits we must take to achieve the "T — > oo" limit eo/h, e' /h —* are vreo, vre' — > (AE). The 
relation between the eQ/h,e'/h — * limit and the T^s oo limit is explicit in a formula 
proven by Kay |]33| : 



/2 fT /2 
^^'g-e'lt'l/n. = lijn - / rft' ■ . (65) 

T^oo T J-T/2 

We consider the correspondence limit of Eq. (0) for the case of n 7^ m, as well as for 
n = m. Consider first n ^ m. Performing the limits as outlined above and making use of 
Eq. (|65D we obtain 



15 



)(n,m|xo)= lim ^5,,{E,,^^ - H{^,)) lim ^ / c?t'e-^'l*'l/V"'-*'5(xo - X(x, -t')) 
= lim ^5,„(E„,„-i7(xo)) lim ^ rftV"— *'5(xo-X(x,-^^^ 
= lim ^5,„(E„,„-/7(xo))2^ rrftV^"-*'5(xo-X(x,-t')), 



1 /""^ 

= dEdX5{En,m- H{{^q))— rftV^"''"*5(xo - X(x, -t') (66) 

Here we have interpreted the limit of 1/T to be dXj^-R and vreo = (Ai?) (ii^^. To see 
that this is correct recall that we have taken the limit as vre' {^E) which corresponds to 
letting T T^s- The inverse of the density of states time can be roughly interpreted as the 
average nearest neighbor Liouville frequency divided by 27r. The average nearest neighbor 
Liouville frequency {AE)/h can evidently be interpreted as dX. Given Eq. (133) we have 



lim (x|ra,m)(n,m|xo) = dEdX T^^ a„.™(x;xo), (67) 

hence proving correspondence. 

Note that in establishing Eq. (^) we have also shown that individual quantum Liouville 
eigenstates (x|n, m), with n 7^ m, do not have correspondence limits when the classical 
system is chaotic. We may also infer the reason: individual states (x|n, m), n ^ m, possess 
essential singularities which cancel in the product (x|n, m)(n, m|xo) to give a well defined 
correspondence limit. 

Consider now the case where = m in Eq. (H). Using Eq. (H) and the fact that 
the classical dynamics is ergodic, we can replace the time average by a phase average and 
rewrite the integral in Eq. (|6^) as 



dt'e-''\''\'^ 5(xo - X(x, -t')) = lim - / dt' 5(xo - X(x, -t')) 

T^oo T J-T/2 

^ 5{E~H{^)) 



(68) 



/rfx'5(E-/7(x'))' 

Substituting this expression into Eq. ( p^ and again interpreting vreo ~ (A-E) ~ dE we 
obtain [Eq. (PI)], the desired correspondence: 
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lim (x|n, n|xo) = (ii? Ts^(x; xo). (69) 

Since 

T^(x;xo) = p:^(xo)p,,(x) = (x|E)(E|xo), (70) 

Eq. ([70|) implies that a product of stationary quantum Liouville eigenfunctions goes to a 
product of stationary classical Liouville eigenfunctions. However, this is not the case for the 
nonstationary projectors. That is, 

T^,a(x;xo) = Ept,A(xo)plA(x) = E (x|i?, A,£)(E, A,£|xo) (71) 

e I 

is not a simple product of Liouville eigenfunctions, but rather a sum of products. Thus Eq. 
(^) implies that the correspondence limit of a product of nonstationary quantum Liouville 
eigenfunctions is a sum of products of nonstationary classical Liouville eigenfunctions, due 
to the degeneracy of the classical states. 

B. Correspondence for the Liouville Spectrum 

In addition to the limit relations for the spectral projection operators discussed in Sec. 
Ill A| , correspondence requires that the quantum spectrum reduce to its classical analog in 
the /i — i> limit. Since the classical Liouville spectrum is continuous we examine spectral 
densities, rather than individual Liouville eigenvalues. The quantum Liouville operator 
5(A — L) can be expanded on the Liouville eigenbasis as: 

5{\ — L) = E 1^5 m){n, m\6{X — L)\n, m){n, m\. (72) 

n,m 

The trace -D(A) of 6{\ — L), is the quantum Liouville spectral density, i.e., 
D{X) = Ti[6{X - L)] = E {n, m\6{X - L)\n, m) 

n,m 

= E {n,m\5{\- \n,m)\n,m) 

n,m 

= Y.^i^-Km). (73) 

n,m 
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With a view toward investigating the classical limit we note that we can rewrite the first 
equality of Eq. (|73D, by inserting the identity 6{X — L) = (27r)"^ / dt e*^'^"^)*, as 

1 pea 

D{\) = Y.7r dt e'^' (n, m\e-'^'\n, m). (74) 

n,m "'-oo 

But, inserting the closure relation ([T8|), and noting that I/|x) = L(x)|x), gives 

(n, m|e"'^*|n,m) = j diL (?T,,m|x) e~*'^^''^*(x|?7,, m). (75) 
It follows that Eq. (^) can be rewritten in the form 

D{\) = Y. — r d^ / d^ in,m\^) e-^^W*( 



X n, m 



n,m 



(n, m|x)(xo|n, m) 



e-'^W*(5(x - xo) 



1 /''^ /" 

= — / dt e'^* / rfxc/xo 
27r j-oo J 

= — dt e'^' / dxrfxo5(x - xo)e-*^W*5(x - xq), (76) 

where we have used J2m,n '"^|x)(xo|n, m) = 5(x— xq). Formally expanding L{x.) in powers 
of Planck's constant and taking the h limit, gives L{x.) — > Lc{x.), where Lc{x.) = iH{x)a 
is the classical Liouville operator. It follows that as /i — >^ 

DiX)^ — / dt e'^' / rfxdxo 5(x - xo)e-*^^('')*(5(x - xq) 

1 f 

^ — dt e^^* / dxdxo S{x - xo)(5(X(x, -t) - xq) 
27r J-oo J 

— / dt e'^' / c/x (5(X(x, -t) - x). (77) 

27r j-oo J 

But 

— / dt e^^* / c/x 5(X(x, -t) - x) = / c/E / rfAo 5(A - Aq) / c/x Ts.a(x; x) 

27r J-oo J Jo J-oo J 

= D,i\) (78) 

where -Dc(A) is the classical Liouville spectral density p5| . 

Thus, in a formal sense we have correspondence, i.e., -D(A) —>■ Dc{X), in the h —>■ limit. 
However, the proof is unsatisfactory because the classical spectrum is highly degenerate and 
this limit is not well defined, i.e., -Dc(A) = oo. This arises from the fact that the classical 



Liouville spectrum is infinitely degenerate due to its stability with respect to variations with 
energy p5|j3^ . In addition, this formal proof provides little insight into the way that the 
spectra approach one another. 

This problem can be bypassed by considering the Liouville spectral density for energies in 
a classically small, but quantum mechanically large, energy interval EQ — t/2 < E < EQ + e/2. 
That is, we define 



D,iEo;\)= E 5i\-\n,m) 

fEo+e/2 

= h dE d{E + hX/2)d{E -hX/2), (79) 

JEo-e/2 

where d{E) = J2n^{E — En). Then, as shown below, expression ( [79|) has a well defined 
classical limit, i.e., 



limlimD,(Eo;A) = D,(Eo;A) (80) 

with 

1 r 

D,{Eo;\) = Si\) + — dte'^' c/x 5(X(x, -t) - x). (81) 

ZVr J-oo JH{yL)=Eo 

Here Dc{Eq; A) is the classical Liouville density of states piJSB] on the energy surface H{x.) = 



Eq. Note that the factor 5(X(x, —t) — x) in Eq. (^) is nonzero only for points x which lie 
on periodic orbits of period t. Thus the integral Jh{x)=Eo '^(^(x, —t) — x) can be written 
as a sum over periodic orbits |34], giving |25,3^ 



j,(E„;A)=^(A)+iE rif:'r'f;™l . m 

Tc ^ kj\ det {Mj{Eo) - I)\ 

where Tj is the period of periodic orbit j, kj is its winding number, Mj is its 2s — 2 x 2s — 2 
dimensional stability matrix, and the sum is over positive traversals of the periodic orbits. 

To show Eq. (|80|) we employ Gutzwiller's formula d{E) = d{E) +dosc{E) for the density 
of states, in Eq. ([791). Here d{E) = {AE)~^ is the average density of states, and dosc{E) is 
an oscillatory correction given by the formula 
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t/osc(i?)~E4c(^) (83) 

j 

where 

. ZM"-(^iW/'' + %). (84) 
y'|det(Mi(E)-/)| 

Here Sj{E) is the action of the periodic orbit j and 7j = o"j7r/2, where cTj is the Maslov 
index of the orbit |36|^8|. The sum in Eq. (|83D is over positive traversals of the periodic 
orbits. 

The Gutzwiller formula for the density of states is not generally convergent, but can be 
made so by broadening over energy ||3^. That is, we replace d{E) by the energy broadened 
density df^^E) = J2n^^l{E ~ ^n) where VL^{x) = 1/fi for — /i/2 < x < /i/2 and is zero 
otherwise. Note that lim^^o ^/^(a;) = S{x). The energy broadening modifies the standard 
Gutzwiller expansion by damping out contributions from very long periodic orbits. 

We rewrite Eq. (^) in the form 

rEo+e/2 

D,{Eo;\)= lim h dE dJE + hX/2)dJE - hX/2) (85) 

^l^(AE) JEo-e/2 



in order to employ the energy broadened (and hence convergent form of Gutzwiller's 
formula. The correspondence limit is now h ^ followed by e, — > (AE) |^ . 

We now separate the Liouville density D^{Eo] A) into its diagonal and off-diagonal parts. 
Since the system is chaotic we assume that it exhibits level repulsion, i.e., that the probability 
of two neighboring energy levels exhibiting an accidental degeneracy is zero. It therefore 
follows that the only contribution to the Liouville spectrum at A = is from the diagonal 
(n = m) terms, 

fEo+e/2 

limA(£'o;A)~ lim h dE Y^VlJE + nX/2 - En)nJE - hX/2 - E^), (86) 

A^O ^j.^{AE) JEo-e/2 „ 

SO that we can rewrite Eq. ( P5| ) in the limit h ^ (i.e., hX ^ for all A) in the form 

pEo+e/2 

DJEo]X)r^ lim h dE ynJE + hX/2- EjnjE -hX/2- En) 

tM-*{AE) JEo-e/2 n 
/■Eo+e/2 

+ Mm h dE [dJE + hX/2)dJE-hX/2)-dl{E)]. (87) 

Ai^(A_B) JEo-e/2 ^ 



Ix^(AE) JEo-e/2 
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Note that as fi {'^E) 

^,.{E + nX/2 - En)n^{E - h\/2 - En) ^ biUX) d{E) (88) 



and so the first term in Eq. (|87|) becomes 

/■Eo+e/2 

lim n dE y njE + hX/2- E„)nJE - hX/2 - E„) 

tJ.^{AE) JEo-e/2 „ 
1 rEo+e/2 

~ he5{h\)- / dE d{E). (89) 

e JEo-e/2 

If we choose e ~ /i/Tmin, where Tmin is the period of the shortest periodic orbit of energy 
Eq, then 

d(Eo) - - / dE d{E) (90) 
and so we may write Eq. (^7[) in the form 

fEo+e/2 

D,{Eo] A) = he6{hX)d{Eo) + lim h dE [d^{E + hX/2)d^{E - hX/2) - dl{E)]. 

fi^{AE) JEo-e/2 

(91) 

Now we assume that d^{E) ~ d{E) + dosc,fj.{E), with ^ JEo^e/2 ^E dosc,fi{E) ~ 0, and substi- 
tuting this expression into Eq. (pT]) we obtain 



DJEq; X) ~ he6(hX)d(Eo) + hed ^(Eq) ~ h lim / dl(E) 

fi^(AE) JEo-e/2 ^ 

/■Eo+e/2 

+ lim h dE doscAE + fi^/2)doscAE-^^/2). (92) 

fj.^{AE) JEo-e/2 

The form of the energy broadened density is now introduced. In particular , 

~ d{E) + yA,,{E)cos{S,{E)/n + ^,), (93) 

where 

T^m^i^^mm^ (94) 
A;,7r;i^|det(Mj(E) -/)| 
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and sinc(x) = sin(a;)/x is the damping function. 

With dosc,fi{E) ~ J2j cos{Sj{E)/h + jj) the last term of Eq. ( p2D can be written 

as 



/ dE d,,,^^{E + h\/2)d,,,^^{E - hX/2) 

JEo-e/2 

J2 / dE A^j{E + hX/2)A^j,{E - hX/2) 

3,3' -^^0-^/2 



cos{Sj{E + nX/2)/n + 7j) cos{Sj>{E - hX/2)/h + 7jv) 

rEo+e/2 



~ ^ / dE A^j{E + hX/2)A^j,{E - hX/2) ■ 

3,3' -^^0-^/2 

{cos [{S,{E + hX/2) - SjiE - nX/2))/h + - -fj,] 
+ cos [{Sj{E + hX/2) + Sj>{E - nX/2))/h + + 7^-/]} (95) 

Now let h 0, followed by e,/i (AE). Observing that - En)^t,{E - Em) ~ 

VL^{E — En)Sn,m/ fJ' if is Sufficiently small (a consequence of level repulsion), if follows that 
d^{EY ~ d^{E)/fi. Taking /i (AE') we see that 



/ rfE rf2(E) ~ / rfE d{E)/{AE) ~ erf ^(^o)- (96) 

JEo-e/2 JEo-t/2 

Substituting this result back into Eq. (^) we obtain the result that 

D,(Eo] X) ~ heS{hX)d{Eo) + lim h dE dose JE + nX/2)doscJE - hX/2). (97) 

fi^{AE) JEo-e/2 



It remains to evaluate Eq. (p5|), a complicated procedure outlined in Appendix B. Then 
taking the limit as e {^E), and substituting the contributions from Eqs. ( |142| ) , ( |144D , 



( |148| ) , and ( |154| ) into Eq. ( pSj) and substituting Eq. ( pS] ) back into Eq. (0) gives the result: 



1 ^ Z7 ^ ^^^t Sj{Eo+\n/2)-Sj{Eo-hx/2) ^ 

D.iE,; X) r. 5iX) + i E ^ ' ' 



ttY fc,|det(M,(Eo) -/)| 



+ ©(/i^-^e^^/'^). (98) 
For s > 2, we thus see that 

... >i. A(^o; A) ^ .(A) . i E ggl-g;*!;" ^ ^=(^0: A). (99, 
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That is, the quantum Liouville spectrum properly approaches the classical Liouville spec- 
trum as h ^ 0. Note that Eq. (pQ]) emerges from Eq. (pSf ) via elimination of essential 
singularities. 

IV. SCAR CORRECTIONS 

In the last section we began with Eq. (p2|) and utilized Berry's formula [Eq. (|46|) ] for 



W{'x;E,e), neglected scar corrections and arrived at a proof of Eqs. @ and @, i.e., the 
correspondence rules discussed above. We now consider the corrections to these limits which 
arise due to the scars from the periodic orbits. 

Consider first the following formula for W^^^^j.^'X] E, e), the scar contribution to iy(x; E, e) 
from periodic orbit j: 

1^4, (x; E, e) ~ =e-^^^/^ ■ 

V|det(M,+/)| 

cos{[S, - ^[JiM, - I)/{M, + I)]l]/h + 7,} ■ 

dt e-l*l/^e^{(^-^^)*-M-Axtn. (100) 

-oo 

Here the variables $, are the 2(s — 1) coordinates of the surface of section transverse to the 
periodic orbit, andxAx= |VF(q)p/m + (p • V)^V"(q)/m^. The derivation of Eq. (|10CI|) is 
given by Berry although he does not write it out explicitly. 

Let /i ^ with e ^ /i/Tmin- [For convenience we drop the j subscripts on Mj, S'j, 7j and 
Tj.] Note that the term in the time integral of Eq. ( |100|) can be neglected because only 
short times count for small h. Thus 

WLA^-. E, e) ~ , '^[^ e-^^'^ ■ 
|det(M + /)| 



cos{[S{E) - i[J{M - I)/{M + m]/h + 7}5,(E - if(x)) 
^|det(M + J)| 

^g-i{[5(£;)-|[j(M-/)/{M+/)]|]/;i+7}]/2 (101) 
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Consider now factors like e^''*^''^'^^ +!)]$, /Tt- the /i — > limit. Note that the integral 

j dz e±'"^^/"/(z) (102) 

[where z = (^i, . . . ,Z]s[) and f2 is an x dimensional matrix independent of z] in the 
limit a can be evaluated by stationary phase to give 

/ dz e±'"^^/"/(z) ~ J2^£L/(0)e±'"^*5"^/^ (103) 
■/ldet(fi)| 



This suggests the existence of a distributional identity 



|det(fi)| ^. 
avrr / 



This formula, when applied to the exponent in Eq. (|101| ), gives 



^'"^^ ^ (106) 



Ydet(M- J) 

[here a = 7rsgn[J(M - /)/(M + /)]/4]. Substituting into Eq. (|10TD gives 

W^i_(x; 6) ~ , ^ ^ e-^/^ cos{5(i?)/;^ + 7}5e(i5 - if(x))5(^(x)), (106) 

^|det(M- J)| 



which is the classical limit obtained by Berry 

To obtain the scar corrections to the Liouville spectral projectors we insert Eq. (^ 



into Eq. (|52D and use Eq. (|106|) for the scar term. Their are two types of corrections of 
the forms: ^^6^{E — H)x scar" and "scar x scar". The latter are of much higher order in h 
and are neglected. The correction terms to Eq. ( ^4]) due to periodic orbits of period T and 
energy En are then 

5„,„(x; xo) = 27r^6i62/i- ' ^_,^T/H ^^,{siEn)/h + 7} ■ 



|det(M-/)| 

J rfyrftirft2e-'^l*^l/^e-^2l*2l/^e2"^(^-^«) -^-y5(^(x + /iy/2)) ■ 

exp{^(E„ - if(x + hy/2))ti/h} exp{t{Em - H{^o - hy/2))t2/n}. (107) 

We neglect h corrections to ^, i.e., we assume that ^(x + hy /2) ~ ^(x). This can be 
justified as follows: (a) expanding ^(x + hy/2) to first order in h and using the fact that 
J(M — I)/{M + I) is symmetric [0, and Eq. (|105|) , allows us to show that 
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5(C(x + hy/2)) ~ 5(^(x)) cos{27r[y • ^] ■ [J(M - /)/(M + /)] ■ (108) 



and (b) noting that the right hand side of Eq. ( p,08| ) is zero unless ^(x) ~ 0, and that the 
argument of the cosine factor is proportional to ^(x) allows us to replace the cosine factor 
by unity. 

Using the short time expansions of Eqs. (|58D and (|59D and doing the integrals over u 
and V gives: 



=e-^^^/'*cos{5(E„)A + 7}5(^(x)) 



|det(M-/)| 
j rftic/t2e"'^l*^l/^e-^^l*^l/^exp{i(E„, - H{^))ti/h} ■ 

exp{t{E^ - i7(xo))t2/^}5(X(x, -ti/2) - X'(xo, -^2/2)). 



(109) 



Changing variables to t' = {ti — t2)/2 and to = and defining eo = ei/2 and e' = 2e2 as in 
Sec. 1111 A| we obtain: 



s-l 



<Sn,m{'^', Xo) ~ ■ , 

2^|det(M-/)| 
(5(^(x)) lim e'/h f dt'e-''\''\/^exp{iXn,mt'}S{xo-M^,-t')). 

For the stationary case n = m we interchange limits via Eq. (^) to obtain: 

5„,„(x;xo) = e-^°'^/^cos{^(E„)A + 7} • 

hJ\det{M-I)\ 



5{En - H{x))5{En - /J(xo))(5(^(x)). 



(110) 



(111) 



The other cross term, due to periodic orbits of energy Em-, gives a similar contribution 
'5„,„(xo; x), and so the overall correction to Eq. (|5^ ) with n = m, due to the periodic orbits, 
is (see also Berry [Q) 

\2 



5„,„(x; Xo) + 5„,„(xo; x) ~ ^ , ^^^^ =e-'oT,/n cos{S,{E^)/h + 7,-} ■ 

, hJ\detiM, - I)\ 



5{E^ - Hi^oMEr, - H{x)) [5(^,(x)) + 5(^,(xo)) 



(112) 



This result may be rewritten in terms of the distributions 



25 



n,o(x; xo) = ^5{E - Hi^oME - if(x))5(^^.(x)). (113) 

which we defined in Paper 1 [^. and which are stationary spectral projectors with uniform 
density on the periodic orbits. Using Eq. (|113|) , Eq. ( |112|) becomes 



T,{AE) 



2 



5„,„(x; Xo) + 5„,„(xo; x) = ^ ^i\^^^_^-eoT,/n cos{Sj{En)/h + 7,'} 

j kjhJ\det{M, - 1)\ 



n,^,o(x;xo) + Ti,^,o(xo;x)J. (114) 

This result shows that the scar corrections to the limit as /i — > of the spectral projec- 
tors corresponding to stationary states are comprised of weighted sums over the stationary 
classical projectors which have uniform density on the classical periodic orbits. For the non- 
stationary case n m, again interchanging limits via Eq. ( |65D and considering only points 
on the periodic orbit, it follows that the integral in Eq. (|110|) can be written as: 



lim e'/h I dt'e-^'l*'l/^e'^-'"*'5(xo - X(x, -t')) 
1 sin{(2fc + l)A„,™r/2} 



= lim 

fc^oo 2k + 1 sin(An^mT/2) 
- citV^"''"*'(5(xo - X(x, -t')). (115) 

r J-r/2 

[Here r = Tj/kj.] Since it can be readily shown that 

lim ^ sin{(2A; + l)A„,^r/2} _ I 1 if Km = 27r//r for / G Z 
fc™ 2k + 1 sin(A,,^r/2) " | q otherwise 

it follows that there are scar corrections only for the nonstationary distributions whose fre- 
quency matches an integer multiple of the frequency of one of the periodic orbits. When this 
condition is met the scar correction to the nonstationary spectral projector (n, m|x)(xo|n, m) 
is obtained from Eqs. (|110|) and ( |115| ) as 



5„,„(x;xo) ~ i^E)h^ cos{S{En)/n + 7}4o(^n,m - g(x)) ■ 

V'|det(M-J)| 

5(^(x))- rftV^"-'"*'5(xo - X(x, -t')). (116) 

T J-t/2 



r/2 
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Similar corrections iSm,n(xo;x) arise from periodic orbits of energy Em, i.e., iS„,m(x;xo) is 
given by Eg. ( |116| ) with En replaced by Em- 

Consider now these nonstationary corrections in more detail. The scar contribution [Eq. 



T6| )] involves the following product of factors 

5,„(K,™-//(x))(5(^(x)). (117) 

The variables ^(x) are effectively zero on a local family of periodic orbits with energies close 
to En and periods close to T. The distribution 5(^(x)) is thus zero except on this local 
family. If the energy En,m lies outside of the neighborhood in energy of the local family 
then the product of delta functions [Eq. (|117| )] will be everywhere zero. As a consequence 



the product is generally zero for n ^ m. The same considerations hold for the scar term 
S'„,m(x; xo) with periodic orbits of energy Em- Thus, scar corrections to the nonstationary 
Liouville eigenfunctions are typically negligible in the semiclassical limit. Only periodic 
orbits of period Tj/kj = 27r//A„^m7 ^ ^ Z, contribute and of these only the ones with energy 
En or Em close to En^m make a nonzero contribution. 

Thus, we see that the stationary and nonstationary contributions of periodic orbits [Eqs. 
( p.l4| ) and ( |116D ] at most make corrections of order fi'^^-^^-^T/n ^Y^^da. vanish in the corre- 
spondence limit (/i — * followed by vre (AE)). Furthermore, in the classical limit these 
corrections are only supported on the measure zero set of periodic orbits. 

V. CORRESPONDENCE: APPLICATIONS TO MATRIX ELEMENTS 

The results obtained above allow us to systematize and extend previous results on the 
classical limiting forms of matrix elements of quantum observables. With our normalization 
of the quantum Liouville eigenfunctions, matrix elements satisfy the following relationship 

{n\A\m) = h-^/' I dx p-„(x) A-(x). (118) 

As a consequence we may relate matrix elements of observables to the spectral projection 
operators via the expression: 
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\{n\A\m)\^ = y^xcixo A'"*(x)(x|n,m)(n,m|xo)A"'(xo). (119) 

We focus on the correspondence limit of matrix elements for observables A with a well 
defined classical analog, i.e., for A^{x) —>■ A{x.) where A{x.) is the classical analog. 
Consider first that Eqs. (|rT|), (|3|), and (gD imply that {n ^ m) 



\im\{n\A\m)\' = '—-^ / dt'e^'-'-' (A(0)A(t'))i.„,. (120) 

h^O h J-oo 

where (^(0)A(t'))£;„^ denotes the microcanonical average of yl(x, 0)y4(x, t') at an energy 
En^m- Noting that 

rftV^"-*' (A(0)A(t'))E„,„ = 27r(A)|„_,„5(A„,™) + 

dtV^-™*' ([A(0) - {A)E^J[A{t') - {A)e^J)e^,^ (121) 
rftV^"-*' ([A(0) - (A)^„,J[A(t') - (122) 



-oo 
oo 



oo 

oo 



if ^n,m. 7^ 0, it follows that 



\im \{n\A\m)\^ = r rftV^' ([^(0) - (A)OKO " (123) 

h^Q h J-oo 

for the nonstationary case. This result is essentially in agreement with the predictions of 
Feingold and Peres |^3[. Their result differs from Eq. ( |123| ) insofar as they are missing a 



oo 



factor of h ^, and they left the energy En,m unspecified. 

The arguments leading to Eq. (|123[) can be readily generalized to the case of mixed 



operators with the result that 
]im{m\A\n){n\B\m) = T rftV^.,™*' ([^(0) - {A)j,^J[B{t') - (5)^„,J)^„,„. (124) 



This important result establishes a connection between the Liouville spectrum of classical 
time correlation functions and quantum matrix elements for chaotic systems. 

The case of n = m also results from Eq. ( |119| ) and Eq. (^Tj) to give the well known 
results: 



lim|(n|i|n)|2= (A) 



En 



\iMn\A\n){n\B\n) = {A)eAB)e„. (125) 
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Scar corrections to the off-diagonal matrix elements are negligible. Scar corrections for 
the stationary case |(r;,|A|n)p are of the form 

E7-#==— e--^^/'cos{^,(E.)/^ + 7,} {A)e. (A)i,„,o, (126) 
j kjhJ\det{Mj — I)\ 

where 

{AYe^o^- r^'dt' A{X{^=,t')) (127) 

Tj J-r,/2 

is the average of A over the periodic orbit j and x-' is a point on the periodic orbit, while 
corrections to Eq. ( |125| ) are of the form 

E , , e-°^^/^os{^,(E.)A + 7.}[ {A)e. {ByE.,o + {B)E„ (128) 

j kjhJ\det[Mj — I)\ 



Comparing Eq. ( |1 23|) to Eq. ( |126|) it is readily apparent that squared off-diagonal matrix 
elements for a chaotic system are of the same order of magnitude as the fluctuations in 
the squared diagonal matrix elements, i.e., 0{h^~^). This confirms an early conjecture by 
Pechukas 

Finally, we compare the chaotic case to the well known ||^ , classical limit for matrix 
elements of individual matrix elements for integrable systems, i.e., 

lim(n|i|m) = -i- / dO A(I„,^, 0) e*("— )-^. (129) 

The most significant difference is that Eq. ( |129|) , which holds for both the stationary and 
nonstationary integrable cases, gives a classical limit for individual matrix elements whereas 



the chaotic results allow classical limits only for products of matrix elements [Eqs. ( |123D , 



( |124| ) and ( |125[ )]. This is a direct consequence of the fact that pn,m has a classical limit for 
regular systems but not for chaotic systems. 



VI. SUMMARY 



Correspondence for chaotic quantum systems has been considered from the viewpoint 
of distribution dynamics in the Wigner-Weyl representation of quantum mechanics. The 
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connections between quantum and classical dynamics have been clarified through the for- 
mulation of the correspondence problem in terms of Liouville spectral projection operators. 
Our demonstration of correspondence for these objects and for the Liouville spectrum shows 
that quantum dynamics is capable of reproducing chaotic classical dynamics in the h ^ 
limit. The mechanism of correspondence here, as in our studies of chaotic mappings, appears 
to be the elimination of essential singularities. Corrections arising from periodic orbits were 
also considered with the result that stationary quantum spectral projection operators have 
corrections in the form of weighted sums of stationary classical projectors on periodic orbits. 
Scar corrections for the nonstationary Liouville spectral projection operators were found to 
be negligible. Applications of our correspondence results to matrix elements revealed con- 
nections between the matrix elements of quantum observables and the spectrum of classical 
time correlation functions. The corrections to the diagonal matrix elements due to periodic 
orbits were shown to be of the same order of magnitude in /i as the square of the off-diagonal 
matrix elements. 



Acknowledgements: We thank the Natural Sciences and Engineering Research Council 
of Canada for support of this work. PB acknowledges stimulating discussions with Professor 
J. Ford, Georgia Institute of Technology. 



30 



APPENDIX A 



We prove identity ( |5DD by showing that the righthand side can be reduced to the left. 
To begin we substitute the definitions of (x|n, n) into the righthand side to obtain 

(q + v/2-q72|n)(n|q + v/2 + q72) • 

(qo - v/2 - q72|m)(m|qo - v/2 + q72). (130) 

Performing the integral over u gives 

h-' J c/vdq'dq" 6{q - qo + q72 - q"/2)e-*(P-P")-^/^e*P 'i'/^e^P" 'i"/'' ■ 
(q + v/2-q72|n)(n|q + v/2 + qV2) ■ 

(qo - v/2 - q72|m)(m|qo - v/2 + q72). (131) 
Now doing the integral over q" gives 

(q- (-v + q')/2|n)(m|q+ (-v + q')/2) ■ 

(qo + (2qo - 2q - v - q')/2|m) (n|qo - (2qo - 2q - v - q')/2). (132) 

Changing variables to v' = — v + q', and v" = 2qo — 2q — v — q' gives a Jacobian determinant 
of absolute value 1/2*. Inserting the definitions of (x|r;,, m) we immediately obtain the 
righthand side of Eq. (^). 

APPENDIX B 

Here we systematically evaluate each of the contributions to Eq. (0), in accord with 
the correspondence limit — > 0, e ^ h/T^i^, followed by e ^ i^E). We assume that Tj{E) 
and Mj{E) vary slowly with energy, i.e., Tj{E±hX/2) ~ Tj{E) and Mj{E±hX/2) ~ Mj{E) 
for small h, and note, as a consequence, that the amplitudes A^^j are also slowly varying, 
i.e., 
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(133) 



and 

A^j,{E + hX/2)^A^^,,{E). 
We thus define the slowly varying amplitudes 

BjAE)^A,,{E)A,,,{E). 
The integrals to be evaluated [Eq. (]95|)] therefore take the form 



Eo+e/2 



dE Bjj,{E) cos 



Sj{E + h\/2) ± Sj,{E - h\/2) 



+ li ± 7i' 



(134) 



(135) 



(136) 



lEo-e/2 

and we must consider three separate cases. 

Case 1: We begin by considering the integrals for which the actions add. Here, because 

we are summing over positive traversals, 
d 



dE 



[Sj{E + hX/2) + Sj,{E - nX/2)] = Tj{E + hX/2) + Ty{E - nX/2) > 



(137) 



and so there are no stationary phase points in the interval Eq — e/2<E<EQ + e/2. It 
can be shown that integrals of the form dz f[z)e^'^^^^/^, with f[z) slowly varying and 
0'(z) 7^ for 2; G [a, 6], can be approximated by 



' dz f{z)e"^^'^/^ 



0'(6) 



0(6)/7i _ 



0'(a) 



J(p{a)/h 



in the limit that h ^ 0. Application of Eq. ( |138| ) to the integral 



Eo+e/2 
Eo-e/2 



dE Bjj,{E) cos 



Sj{E + hX/2) + Sj>{E - hX/2) 

h 



+ 7i + 7i' 



(138) 



(139) 



gives 



rEo+e/2 

/ dE B..,(E) COS 

JEn-e/2 



'Eo-e/2 



Sj{E + hX/2) + S.y{E - hX/2) 



+ 7i + 7i' 



h{ 



B,,,{Eo + e/2) 



sm 



Tj{Eo + e/2)+Tj,{EQ + e/2) 
Sj{Eo + e/2 + hX/2) + Sj>{Eo + e/2 - hX/2) 



n 



+ 7i + 7/ 



B,,iEo-e/2) 



sm 



T^{Eo-e/2)+Ty{Eo-e/2) 

Sj{Eo - e/2 + hX/2) + Sj,{Eo - e/2 - hX/2) 



h 



+ 7j + 7i' 



}• 



(140) 
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Taking the limit as e {AE), Bjj>{Eo ± e/2) ~ Bjj,{Eo) and Tj(Eo ± e/2) ~ Tj^Eo), so 
that we can rewrite Eq. ( |140[ ) in the form 



So+e/2 



2hBjji (^Eq 



dE Bjj>{E) cos 



cos 



SjjE + hX/2) + SjjE - hX/2) 

ti 

SjjEo + hX/2) + SyjEp - hX/2) 

h 



T,iEo)+TAEo: 
sm[{T,{Eo)+T,iEo))e/h] 

\ S,{Eo + hX/2) + SAEo-hX/2) 

~ ZejDj j'ii^o) cos 

[ h 

smc[{T,{Eo)+T,>{Eo))e/h]. 



+ 7i + If 



+ 7i + If 



+ ij + ly 



1411 



Further simphfications show that these contributions take the form 



/ dEB..>(E) COS 

JEn-e/2 



JEo-e/2 

2{AE)Bjj'{Eo)cos 



SjjE + hX/2) + Sj>{E - hX/2) 



+ 7i + 7j' 



SjjEo + hX/2) + SfjEp - hX/2) 



h 



+ 7j + 7i' 



(142) 



Now consider the terms in Eq. (|136|) involving a difference of actions. 

Case 2: Of the terms with j ^ j', some will have stationary phase points in the interval 
Eo — e/2 < E < EQ + e/2, and others will not. For those terms which do not have stationary 
phase points 



'Eo+e/2 

dE Bjj>{E) cos 

So -e/2 

2hBjji (^Eq 



'Sj{E + Xh/2) - Sf{E - hX/2) 



+ 7i - ir 



■ cos 



Sj{EQ + hX/2) - Sj.{Eq - hX/2) 



h 



Tj{Eo)-T,,iEo 
sm[{T^iEo)~TAEo))e/h] 

eos \SAEo + nX/2)-SAEo-nX/2) 



+ 7i - ir 



n 



+ 1] - ir 



sinc[(T,-(Eo)-T,.(Eo))e//i], 



(143) 



as one can easily show by applying Eq. ( |138| ). In the limit as e ^ {AE) we find that 



Eo+e/2 



dE Bjj,{E) cos 



Eo-e/2 

2{AE)Bjj,{Eo)sm 



Sj{E + hX/2) - Sj\E - hX/2) 



+ Ij - If 



Sj{Eo + hX/2) - Sy{Eo - hX/2) 



+ 7i - If 



(144) 
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Consider those terms which do have stationary phase points E' G [£"0 ~ ^/^^ Eq + e/2]. 
We define 0(E) = Sj{E + hX/2) - Sj,{E - hX/2). Expanding (f){E) about the stationary 
phase point E' gives (j){E) ~ + \(t)"{E'){E - E'f. Now consider that 



This approximation resembles the stationary phase approximation except that we have 
retained the original limits on the remaining integral instead of replacing them by ±00. 
Evaluation of the remaining integral in Eq. ( |145| ) then yields |^ 



erf I + ,1, _ E')^ - erf |/^,B„ ^ _ ^0] 



(146) 



In the limit as e — >■ (AE) we obtain 

lEo-e/2 



r^^'^ dE 5,- ,,(E)e^<^(^)/" ~ (AE)fi,- ,,(E')e*^(^')/^*^"(^')(^°-^')'/2^ (147) 



or 



/ dE B., A E) COS 

JEo~e/2 



Sj {E + hX/2) - Sy {E - hX/2) 

h 7j — 7/ 



h 

~ (AE)5j-jv(E')cos($), (148) 

where $ = (j){E')/n + (j)"{E'){Eo - E'f/2h + 7^- - 7^-.. 

Case 3: Finally, we consider the terms in Eq. (|136D where j = j' . First note that there 
are no true stationary phase points, i.e., Tj{E + hX/2) — Tj{E — hX/2) 7^ for A 7^ 0. Second, 
note that the action Sj{E) changes by Sj^Eq + e/2) — Sj^Eq — e/2) = ASj^Eq) ~ Tj{EQ)e 
over the interval Eq — e/2 < E < Eq + e/2. Since e ~ h/T^i^ it follows that ASj{EQ)/h ~ 
27rTj(£'o)/^min- In general Tj{Eo) » T^in and so S/h changes dramatically over the interval 
Eq- e/2 < E < Eo + e/2. Similarly, A[Sj{Eo + hX/2) - Sj{Eo - hX/2)]/h ~ 27r[Tj{Eo + 
nX/2) - Tj{Eo - hX/2)]/T^i^, and since even \Tj{Eo + hX/2) - Tj{Eo - hX/2)\ » T^nm 
in general, it follows that the phase of the cosine factor will oscillate many times over the 
interval Eq — e/2 < E < Eq + e/2. By contrast, over the interval Eq — hkj/2Tj{Eo) < 
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E < Eo + hkj/2Tj(Eo), A[Sj{Eo + h\/2) - Sj{Eo - n\/2)]/n ~ 27rkj[Tj{Eo + h\/2) - 
Tj{EQ-h\/2)]/Tj{EQ). In general |Tj(Eo + /iA/2) -T,(Eo - /iA/2)|/T,(Eo) « 1 for h small 
and so the phase is effectively stationary over the small interval Eq — hkj/2Tj{EQ) < E < 
Eq + hkj/2Tj{Eo) about Eq, and we will therefore break the integral into three parts: 

JEo-e/2 J^Jy y y 

i'Eo+hkj/2Tj {Eo) i-Eo-hkj/2Tj (Eq) 

[ dE + dE 

JEo-hkj/2Tj{Eo) JEo-e/2 

rEo+e/2 SAE + n\/2) - SAE -n\/2)^ , , 

+ / dE\B.AE)co^{^ 1 ) o\ ^ _ 

JEo+hkj/2Tj(Eo) h 

The first term can be approximated via 

fEo+hk,/2T,{Eo) ,SAE + h\/2)-SAE-h\/2)^ 
/ dE B. AE) cos(^^ 1 I _ 

J Eo-hkj /2Tj{Eo) ' Tl 

kjhBj,,{Eo) ^ Sj{Eo + h\/2)-S,{Eo-hX/2) ^ 

cos(^ — ^ ), (150 

while the other two terms can be evaluated via Eq. ( |138D to give 



2hBj^,{Eo) , S,{Eo + hX/2) - SjjEo - nX/2) ^ _ 

COS! ) * 



and 



Tj{Eo + h\/2) -Tj{Eo-h\/2) ^ Ti 
sin [(e/2 + hkj/2T,{Eo)){T^{Eo + Ti\/2) - T^{Eo - hX/2))/n] , (151) 



Eo+e/2 ,SAE + h\/2) - SAE -h\/2)^ 

dE R- AE) cosM^ 1 I J\ _ 

Eo+hkj/2Tj{Eo) h 

2nBj,j{Eo) Sj{Eo + h\/2) - Sj{Eo - h\/2) , 

COS ( 



Tj{Eo + hX/2)-Tj{Eo-hX/2) ' % ' 

sin [(e/2 - hkj/2Tj{Eo)){Tj{Eo + h\/2) - Tj{Eo - h\/2))/h] . (152) 

Combined, these last two terms contribute 

sinc[e(T,(Eo + hX/2) - Tj{Eo - hX/2))/n] ■ 

cos[27ikj{Tj{Eo + \h/2) - Tj{Eq - Xh/2)) /Tj{Eq))]. (153) 
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Therefore, in the Umit as e — (AE) the total contribution for case 3 is 



/ dE B. J E) COS 

JEn-e/2 



Eo-e/2 h 

hkjBj^j{Eo) .Sj{Eo + h\/2) - S,{Eo - h\/2). 

^ UE,) n ) 

cos[2Tik j{Tj{EQ + Xh/2) - Tj{Eo - \h/2)) /Tj{Eo))]. (154) 



In summary, evaluation of the three types of integrals in Eq. (|95D in the correspondence 
limit, yields the contributions of Eqs. ( |142| ) , ( |144| ) , ( |148| ) and ( |154| ) which combined make a 



total contribution of the form 



hk,B,AE,.) ^^^S,(E„ + Xn/2) - S,(E„ - o(,.-,Wa). ,155, 

Tj{Eo) h 
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